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Using 20 months of CPU time on our special purpose computer "Percola" we 
determined the exponent for the normal conductivity at the threshold of three- 
dimensional site and bond percolation. The extrapolation analysis taking into ac- 
count the first correction to scaling gives a value of t/i/ ~ 2.26±0.04 and a correction 
exponent uj around 1.4. 

1. Introduction 

Percolation is the simplest model for a random medium and is applied to 
describe porous media, gelation, mixtures of electrical conductors, and many other 
disorder dominated situations. It is defined on a lattice by occupying with probabil- 
ity p randomly and independently either the sites or the bonds with conductors. If 
one considers that locally current can only flow between conductors that are nearest 
neighbors, then only above a critical value pc a global current can flow through the 
system. The global electrical conductivity S rises like 

^--{p-PcY (1) 

and the accurate determination of t, or more precisely oit/v where v is the exponent 
of the correlation length, is the aim of the present paper. 

In the '80s, after a heated controversy the critical exponents for the conductivity 
of two-dimensional percolation were determined quite accurately and various 
conjectures could be ruled out. The three-dimensional situation remained less 
clear because of the larger error bars. Only recently new efforts have been made '^'^l 
to obtain more precise values for t/v in three dimensions using considerably larger 



computational means. Gingold and Lobb t^l find tjv = 1.21% ± 0.012 and Batrouni 
et al. [^1 claim tjv = 2.282 ± 0.005. 

Using the specially constructed processor "Percola" , assembled in Saclay by 
J.-M. Normand and M. Hajjar between 1984 and 1987, we did exact calculations on 
narrow bars of random conductors ("strip-bar method"), a technique that has given 
very precise values in two dimensions PI and for the superconducting-conducting 
mixture in three dimensions [^1. For the three-dimensional conductivity problem 
treated in this paper this method a dozen years ago gave t/v = 2.2 ± 0.1, however, 
with much smaller computational investment 1^1. 

In the following section we describe the methods, in the next section we present 
our results and analyze them, and finally we conclude. 

2. Method 

The method wc use is a straightforward generalization to three dimensions of 
the technique of Derrida and Vannimenus 1^°] . We calculate the electrical current 
on a simple cubic lattice of width m and depth n in the horizontal directions, and 
of length L vertically, growing from bottom up. Periodic boundary conditions are 
used in the y-direction, defined in fig. 1. Each bond between nearest neighbors is 
randomly conducting, in bond percolation if the bond is occupied (with probability 
p), in site percolation if both neighbors are occupied (each with probability p), 
otherwise it is insulating. The top and the bottom planes are two equipotentials 
(fixed boundary condition). One defines by /j and Vi the current and potential 
at the end of each line i, i = 1 being the index associated with the top plane. 
The conductance matrix aij, given through Jj = CTjjVj, is updated via 

a'^ j = aij — CTj „crajr/(l -|- (Ta,a''') if a longitudinal bond is added at line a and 
via <t[ j = (Ti_j + {5a. j — Sfj.j) ■ {Sa.i — Sp.i)/r if a transversal bond is added between 
lines a and /? where r is the resistance of the bond. In that way the conductance 
is obtained exactly. A very large L (bar configuration) automatically averages 
over many configurations leaving no ambiguity over whether one should average 
the conductivity, the resistivity, or functions of them. This is an advantage over 
methods that average over many disjoint cubes [3>5,6,io] 



Fig. 1 Bar of length L, width m and depth n. We calculate the conductivity between 
the two black planes. 



Wc will consider, as in [^1, the eases rn = n — 1 for bond pcreolation and m = n for 
site percolation. The conductance of the bar per length is then (Ti^\/L which in the 
limit L — *■ 00 goes to a value it„. The calculation is performed for different linear 
sizes n at the percolation threshold pc = 0.248812 [^^1 for bond percolation and 
Pc = 0.311605 '^^1 for site percolation. The critical exponent t/v is then extracted 
by using the finite size scaling: 

cr„ =an-*/''(l + 6n-'^) (2) 

where a and h are non-universal constants and lo is the (first) universal correction 
to scaling exponent. Since i/, the correlation length exponent, is known rather 
precisely {v = 0.875 ± 0.003) I^^l one can determine t. 

3. Results cind Analysis 

The values obtained with "Percola" using a CPU time of 20 months are given in 
table 1 for the case of bond and of site percolation. Percola ^''^ , a 64-bit floating point 
processor runs at 25 MFlops but is specialized to deal with the algorithm described 
in the previous section for which it runs faster than a Cray-XMP mono-processor. 
The lengths of the bars presented in table 1 exceed those calculated before 1^1 by 
about a factor of 1000. 
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Table 1 Conductance an per length for bond (left) and site (right) percolation for 
different widths n of the bar. We also show the statistical mean square deviation 
A(T„ and the length L of the bar for each case. 



We analyzed the data with various techniques. Finding that t/v which minimizes 
the squared error for fixed uj (as done in rcfs. [2] and [8]) gives t/v = 2.15 ± 0.12 
for the sites and t/u = 2.26 ± 0.06 for the bonds. In both cases u ~ 1.4 gave the 
best fits. Since universality states that these exponents should be the same only 
the intersection of the error bars. i.e. t/v = 2.23 ± 0.06, should be relevant. 

The data for bond percolation are closer to the asymptotic regime as can be seen 
from the error bars and the curvature in a plot log(cr„n*/'') against n~'^. Looking for 



bond percolation at successive slopes (local derivative of log((7) vs log(n)) against 
n~'^ gives the most plausible curve for u around 1.4. In that case t/v = 2.25 ±0.05. 

Keeping t/v fixed and finding the uj which minimizes the error (and number of 
iterations) in the Runge-Kutta procedure of minimizing the mean square deviation, 
gives a dependence t/u against to as shown in fig. 2. For w around 1.4 we find 
t/v = 2.26 ± 0.04 for bond percolation. Unfortunately, as opposed to the case of 
the superconductivity exponent the prefactor of the site percolation data has the 
same sign as for bond percolation. Therefore the curves for site percolation in fig. 2 
do not intersect with those of bond percolation which does not allow us to increase 
the precision by assuming universality and combining the two data sets. 



Fig. 2 Leading exponent t/v against the first correction exponent w for bond perco- 
lation omitting all sizes of width less or equal to n;, for n; = 2 (diamonds), 3 (crosses 
+), 4 (squares) and 5 (crosses x). 

4. Conclusion 

Although we have obtained extremely precise data for the conductance of nar- 
row systems our error bars for the exponents t/v extrapolated to infinite sizes are 
not as small as obtained by other methods '^'^l. This is due to the fact that we 
have seriously taken into account the first corrections to scaling which seems im- 
perative for the case of small systems. As with many simulations of high precision, 



statistical errors arc much smaller than the systematic ones. In contrast to the 
two-dimensional case and the superconducting case in three dimensions '^1 the 
relation between leading and correction exponent is not in opposite sense for site 
percolation as compared to bond percolation giving us a particularly unlucky sit- 
uation. We conclude t/v = 2.26 ± 0.04 consistent with previous values and the 
Alexander-Orbach conjecture W . 



We thank D. Stauffer for help in the evaluation of the data and advice and for 
helpful suggestions about the manuscript. 
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